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Abstract: Motivated by the possibility of inflation in the cosmic landscape, which may be 
approximated by a complicated potential, we study the density perturbations in multi-field 
inflation with a random potential. The random potential causes the inflaton to undergo a 
Brownian-like motion with a drift in the Z}-dimensional field space, allowing entropic per- 
turbation modes to continuously and randomly feed into the adiabatic mode. To quantify 
such an effect, we employ a stochastic approach to evaluate the two-point and three-point 
functions of primordial perturbations. We find that in the weakly random scenario where the 
stochastic scatterings are frequent but mild, the resulting power spectrum resembles that of 
the single field slow-roll case, with up to 2% more red tilt. The strongly random scenario, in 
which the coarse-grained motion of the inflaton is significantly slowed down by the scatter- 
ings, leads to rich phenomenologies. The power spectrum exhibits primordial fluctuations on 
all angular scales. Such features may already be hiding in the error bars of observed CMB 
TT (as well as TE and EE) power spectrum and have been smoothed out by binning of data 
points. With more data coming in the future, we expect these features can be detected or 
falsified. On the other hand the tensor power spectrum itself is free of fluctuations and the 
tensor to scalar ratio is enhanced by the large ratio of the Brownian-like motion speed over 
the drift speed. In addition a large negative running of the power spectral index is possible. 
Non-Gaussianity is generically suppressed by the growth of adiabatic perturbations on super- 
horizon scales, and is negligible in the weakly random scenario. However, non-Gaussianity 
can possibly be enhanced by resonant effects in the strongly random scenario or arise from the 
entropic perturbations during the onset of (p)reheating if the background inflaton trajectory 
exhibits particular properties. The formalism developed in this paper can be applied to a 
wide class of multi-field inflation models including, e.g. the N-flation scenario. 
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1. Introduction 



The study of the inflationary universe is mainly on the slow-roll model with a single inflaton 
field. However, motivated by a variety of reasons, including those from super string theory, 
multi-field inflation has received a lot of attention recently. An important reason to focus 
on the single field case is simplicity. A typical multi-field inflationary scenario is clearly 
very complicated and not predictive. There are special multi-field models which are easy to 
analyze while still predictive: models where the fields do not interact with each other during 
inflation were analyzed first in Ref.|l]]. A special case of this class of models is when all 
inflaton fields have the same potential, such as assisted inflation 0. While in more general 
cases, e.g. N-inflationQ, and similar scenarios f§ |], each inflat on can have its own potential. 

If the inflaton fields couple to each other during inflation and reheating, the problem seems 
rather intractable. Motivated by the realization of inflation in superstring theory, we expect a 
multi-field inflationary scenario. This is particularly true if inflation takes place in the stringy 
cosmic landscape §, where the number of moduli is probably as large as hundreds. Even if not 
all of them participate in inflation, it is likely that a large subset of moduli contribute directly 
or indirectly to the inflationary scenario. Their interaction can be rather complicated (maybe 
strongly interacting as well), generating a cosmic landscape that is random (aperiodic) in 
some directions and quasi-periodic (within certain field range) in others. 

A realistic illustration could be the following potential for moduli p\ and axions </>iJ?J, 

Vrifr, <k) = V(j>j) - ^ cos + cos - + . . . + Ufa, • (1.1) 

/j's are the axion decay constants and a« = Mfe~ i^st ^> where Mj is a natural mass 
(say, the string) scale, and the ith. instanton has action S\ nsV V(pj) contains the potential 
coming from the moduli and includes D-term contributions. Presumably, it contains the 
dominant contribution to the vacuum energy density in Vrfa,(f>i), while U(pi,<j>i) expresses 
the couplings of the moduli and the axions. In the region where Vfa) is relatively flat or 
relatively localized, extensive inflation may take place, and the inflaton is represented by the 
collection of <f>i axions. With very small ctj, the wavefunction is either not trapped (because 
binding in higher dimension typically requires a deeper attractive potential) or tunneling out 
of any classical vacuum site is very rapid. This potential is periodic except for the interaction 
term Ufa, fa), which is perturbative in general, but it can also introduce relatively strong 
interactions at times. Here, Ufa, (pi) plays the role analogous to impurities in a crystal or 
doping in condensed matter physics. 

On a more general ground, the complicated scalar potential could be described using the 
random matrices ||. Generically, the potential probably looks rugged at short distances but 
smoother at coarse-grained level. We shall assume the randomness in the potential comes 
in different scales. Roughly speaking, we have in mind the following scenario: the small 
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scale randomness of the potential causes the inflaton to undergo Brownian-like motion in 
the .D-dimensional field space, allowing the entropic perturbations to feed into the adiabatic 
perturbations during inflation. The large scale behavior of the potential allows the inflaton 
to slowly move down in the landscape and eventually end the inflationary epoch. The ending 
can be brought about by the appearance of tachyons, fast tunneling to a low point in the 
potential, or a quick roll down when the inflaton suddenly reaches a steep slope. 

It is important to point outQ that the number of scalar fields relevant for inflation should 
be bigger than 2. Otherwise, the inflaton wavefunction will be trapped (i.e., localized) in 
the inflaton potential. If so, eternal inflation is unavoidable and inflation does not end. Our 
scenario implicitly assumes there is no eternal inflation and this assumption becomes more 
likely as D, the number of scalar fields participating in inflation, increases beyond 2. Here 
we have in mind a relatively large D (say, dozens to hundreds). 

In slow-roll inflation, the large number N e of e-folds of inflation follows from the flatness of the 
inflaton potential. In general, this requires a fine tuning. Here, the large scale structure of the 
inflaton potential has to be reasonably flat too. However, the large scale flatness constraint 
is substantially relaxed compared to that in the slow-roll case, since the inflaton exhibits 
a random walk like motion and scatterings tend to slow down its overall motion. For this 
reason, we consider multi-field inflation in the landscape to be quite natural. 

Motivated by this picture, it is reasonable to ask if one can make any predictions about such 
an inflationary scenario. Although the analysis is rather general under certain conditions, it 
helps to have a concrete example in mind. As an illustration, we could imagine an extended 
version of brane inflation, where a D3-brane is moving in the cosmic landscape. Its position 
during inflation includes both its position in the moduli space and its position within the 
6-dimensional compactified bulk. Its motion includes classical rolling and scattering, which 
may result in percolation. Quantum fluctuation leads to density perturbation. Presumably, 
quantum diffusion may be included as well. One may consider this scenario as a generalization 



of chain inflation[10], which includes only repeated fast tunnelings, and related scenarios in 
Ref.0. 



Note that, in the absence of quantum fluctuation, no cosmological perturbation is generated 
in the classical but random trajectory of the inflaton. Even for a convoluted random looking 
trajectory, the classical path of the inflaton will be the same at every causal patch of the 
universe during inflation and so no cosmological perturbation is generated in the absence 
of quantum effects. So the cosmological perturbation is entirely sourced by the quantum 
fluctuation. As is well known, in multi-field inflation, the quantum perturbation results in 
the adiabatic mode and the entropic modes of the cosmological perturbation. If the classical 
trajectory is a straight line in the field space, the entropic modes and the adiabatic mode 
remain decoupled, and the density perturbation comes entirely from the adiabatic mode. The 
randomness of the classical background allows for the random mixing of adiabatic and en- 
tropic modes, which is how the randomness of the classical background path enters into the 
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cosmological perturbations, and how the resulting power spectrum incorporates the random- 
ness of the inflaton potential. In general, a random potential will lead to a combination of 
classical percolation and quantum diffusion. However, knowing that the density perturbation 
is very small, i.e. dp/p ~ 1CP 5 , it should be a very good approximation to treat the quantum 
diffusion as a quantum fluctuation around the classical (but complicated) background trajec- 
tory. In this case, the data will never be able to reveal individual bumps. Our analysis in 
this paper is carried out in this approximation. A full analysis of quantum diffusion will be 
interesting but is beyond the scope of this paper. (Of course, one may take the point of view 
that, since 5p/p is so small, such an analysis may not be necessary.) 

Here we would like to show that if the potential is random enough, some definitive statements 
about such a multi- field inflationary scenario are possible. In fact, the predictions can be quite 
specific and detail independent. 

• Even if the inflaton moves relatively fast, its drift velocity can be small, dictated by the 
slope of the potential over large field scales. Inflation naturally lasts for many e-folds. 

• One may easily recover the usual slow-roll, almost scale-invariant power spectrum, with 
a slight additional red tilt (up to 2% when scatterings are frequent but mild) in the 
power spectrum index. 

• The running of the power spectral index can be large, because the conversion of entropic 
modes into the adiabatic mode introduces additional scale dependence. 

• The tensor to scalar ratio r can be large. On one hand, it is suppressed by the amplified 
scalar power spectrum; on the other hand, it gets enhanced by a factor equal to the 
square of the ratio of the Brownian-like motion speed to the drift speed, which we expect 
to be much greater than one in the strongly random scenario. The net effect may still 
be an enhancement in r. 

• Non-Gaussianity is generically suppressed by the growth of the adiabatic mode after 
horizon exit due to the presence of extra light fields. However, it may get enhanced by 
resonant effects E2[. On the other hand, the large number of strongly coupled fields can 
enhanced the non-Gaussianity induced during the onset of (p)reheating. 

• We generically expect random fluctuations in the scalar power spectrum due to scat- 
terings of the inflaton. Such fluctuation may appear with random amplitudes and 
periodicity. The precise shape of the power spectrum is not quite informative here, 
but the variance of such fluctuations is well determined by the underlying microscopic 
properties of the scalar potential (landscape). The currently observed TT power spec- 
trum, looking at the resolution of each multiple moments, is not smooth. The error bar 
of each data points (roughly 10%) is comparable to the fluctuation of data itself (for 
100 ~ ^ IS 800) so it is hard to claim any fluctuating features in the currently ob- 
served power spectrum. Furthermore, current analysis based on binned data points from 
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a few tens of multiple moments may have smoothed out such fluctuations. However, 
with more observational data coming in the future, we expect that analyzing the power 
spectrum at the resolution of single multiple moments is possible and the fluctuations 
in the power spectrum can certainly be detected or falsified if the error bars become 
substantially smaller than the primordial fluctuations. Furthermore, if such fluctuations 
are detected in the TT power spectrum, one expects to see the same pattern in TE and 
EE correlations as well. 



The paper is organized as follows. The relevant properties of the multi-field inflationary 
scenario are reviewed in Section ||. In Section |3|, we consider the model where the randomness 
in the potential leads to a Brownian-like motion for the inflaton. That is, only classical 
percolation is included. The power spectrum is discussed in Section || and the non-Gaussianity 
is discussed in Section |5[ In Section ||, we discuss the more general situation where the random 
scattering is also caused by quantum diffusion. We attempt to argue that the generic results 
obtained in Section ^| and Section |5| may hold even in the presence of such quantum diffusion. 
Section ^ contains the summary and some remarks. Some details are relegated to the two 
appendices. 



2. Review of the Multifield Inflation Formalism 



In this section, we review the multi-field inflation models for D light canonical scalar fields. 
We start with the action 




where we define X = — \ g^ v d^cj) 1 d v cf)i, with I = 1, 2, . . . , D. The 4 dimensional Friedmann- 
Lemaitre-Robertson- Walker metric takes the form 



ds 2 = g iiv dx il dx v = —dt 2 + a(t) 2 dxi dx l . 



We have set JWj = 8vrG = 1. 

In this paper, we will mainly consider the canonical kinetic term for simplicity. Generic 
kinetic terms could be a function of X = — ^Gjj(4>)g flu d ll (j) I du4 >J \l^\: or more generally 
(especially for the motion of D-branes) as shown in Ref. J15||, the full action should be a 
function of X IJ = — ^ g^ v d^cj) 1 d v c/) J , and involves non-trivial terms like XjjX ij . To simplify 
the analysis, we consider only canonical kinetic terms. The formalism developed in this paper 
could be generalized to the scenarios with non-canonical kinetic terms. 
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2.1 The Background Equations 



The background equation of motion for the scalar field (j) 1 i 



is 



4>i + w<t>i + Vi = o 



ay 



(2.2) 



0^ ' 



We shall consider both the short distance scale properties and the large distance scale prop- 
erties of this multi-field inflaton potential V ((/)). The coarse-grained potential will be denoted 
by V, which leads to the drift velocity of the inflaton. The scattering by the short distance 
property of the potential will lead to Brownian-like motion. 

The background Friedmann equations are 



In the presence of multi-fields, it is useful to define a set of orthonormal basis {e n } in field 
space. The first vector is defined to be along the direction of the background trajectory, 



The rest D — 1 vectors are orthogonal to e a and they are collectively called e s . Among all 
the e s , it will prove useful later to choose one along the direction e^. This particular unit 
vector is denoted by e K , and it points along the curvature radius of the background path. 

The background equation of motion can be projected to the e a direction, and we get 



One should note that the basis {e n } defined here does not point to fixed directions in the 
field space. They are defined on the background path of the inflaton, and their orientations 
are time dependent. 

2.2 Linear Perturbations 

In the longitudinal gauge, the metric perturbation can be written as 



In this particular gauge, the gauge invariant curvature perturbation £ and scalar field pertur- 
bation Qi takes the form, 





& 2 = 2X . 



a + 3H& + V a = , V a = Vi^ a . 



(2.3) 



ds 2 = -(1 - 2<$>)dt 2 + a(t) 2 (l + 2^>)dxi dx { . 



(2.4) 
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C and Qi are related through the Einstein's equation, 



On super-horizon scales, the right hand side is negligible, and we have 



C ^ --Q* , Qa = eiQ! . (2.6) 

(7 



According to Ref.[16], £ evolves with time on super- horizon scales, 



H p 

C = 1 — <5Pnad , Sp nad = 5p- -Sp . (2.7) 

P + P P 

where Sp na d is the non-adiabatic pressure perturbation. 

For single field inflation, 5p na( j = 0, so ( is conserved. However, in the case with D canonical 
scalar fields, on super-horizon scales where the spatial gradient term can be ignored, Ref.[17] 
used a hydrodynamic approach and found that 

C = —eiQj (2.8) 
a 

One sees that ( is generically not conserved when the background path makes a turn, i.e. 
e ff 7^ 0, and it is only sourced by entropic perturbations along the e CT direction. 

2.3 The SN Formalism 

The primordial curvature perturbation £(£, x) generated during multi-field inflation can also 
be described using the SN formalism^], which not only reproduces the result in the previous 
section at linear order, but naturally allows us to go to second order perturbations as well. 
The main idea of the SN formalism is that during inflation, separate Hubble volume evolves 
independently, and they differ by how much they have expanded relative to each other towards 
the uniform energy density slice at the end of inflation. The difference 5N is caused by the 
super-horizon quantum fluctuations of the inflaton field and is related to the large scale 
curvature perturbation £ on the uniform energy density slice through 

( = 5N = N(<t>}(x), 0f (#)) - N[ , (2.9) 

where Nf is the number of e-folds on the spatially flat slice, and A^(0|(x), (pf ((p^)) is the 
number of e-folds on the uniform energy density slice. ^J(x) denotes the field configurations at 
the time t* of horizon crossing. cj)f denotes the field configurations at the end of inflation, and 
they generically depend on the initial in the multi-field case. For the slow-roll single field 
case, we have 5N = H5t, which was used to evaluate the scalar density perturbation Jl9|, |20[| . 
More generally, we can expand SN as 

SN = N I Q I + l -Nu {Q I Q J - (Q 7 Q J )) + . . . , (2.10) 
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where we have defined Nj = dN/dcf) 1 , and Njj = d 2 N/d(f> I d(f) J . Q is the perturbation of 
(f) 1 in the spatially flat gauge. For future reference, the perturbation along the direction 
is denoted by Q a and perturbations along the e s directions are denoted by Q s . We have 
N a = iV 7 e^ and N s = iV/ef. 



For the single field case, the leading order term in the 5N expansion gives 

Q ■ (2.11) 



t* 



From the 5N perspective, since different Hubble volumes all follow the same trajectory, the 
only source of 5N is the initial value of (f> they start with. 5N does not change once 5(f) get 
frozen on super-horizon scales, leading to the conservation of £. 

The non-conservation of £ in the multi-field case can also be understood in the 5N formalism. 
We start by writing the e-folds along the unperturbed path as 

N = — da , (2.12) 

Jo* v 

from which we perform the first order expansion, 



5N = -- 
a 



H 



(dja^Q 1 + - {d I a E )Q I 



a 



The first term above comes from the shift in initial value a*, and it corresponds to the 
adiabatic perturbation in the single field case. The second term containing <9/<te arises when 
the hyper-surface of the end of inflation is not orthogonal to the background inflaton path. 
The last two terms are non-local and depend on the complete inflaton trajectory after t*. 
They capture the fact that under entropic perturbations, the inflaton follows a new trajectory 
nearby, and the new trajectory has different length (the third term) and also different speed 
(the fourth term), both contributing to 5N. The last term only contains Q s perturbations, 
as it represents perturbations orthogonal to the integral path. 

The d\OE term is sub-leading when the inflaton executes Brownian-like motion, so we will 



ignore this term for a moment. In Section 5.2, when we look at the non-Gaussianity at the 



onset of (p)reheating, this term will be important. The non-local terms in Eq.( |2.13| ) can be 
evaluated using a geometric approach (see Appendix [A] for detailed derivations), and we get 



5N = -* 
a 



rtE off 

Q a - dt—eiQL (2.14) 
Jt* ° 



The Q a perturbation has the same interpretation as in the single field case, it fluctuates 
tangent to the inflaton trajectory, leading to 5N between Hubble volumes that follow the 
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same trajectory. However, there also exists fluctuations along the direction e CT (perpendicular 
to e CT ), which kick the Hubble volume off the original expansion trajectory. 

The off-trajectory perturbations are exactly the source of the non-conservation of £. In fact, 



using the relation £ = 5N together with Eq.( 2.14 ), we get 



C = -—eiQ I 
a 

which is the same as what one gets from the hydrodynamic analysis Eq.fl2.£|). 

By comparing the 5N expansion Eq. fl2.14Q with Eq.( 2.10]) at first order, we identify 



N a Q a = -—Q a (2.15) 
a 

N s Q s \ t * = -l dt — eiQj. (2.16) 



t* 



a 



It is also useful to write in components, 

tt rtfi 27^ 

N a \ t * = --, N s \ t * = - dt — (e a -e* s ). (2.17) 
cr J t * a 

We see that N s depends non- locally on the whole inflaton trajectory from t* to making 
it difficult to evaluate if the background inflaton trajectory makes several turns. 

We should also point out that unlike the adiabatic perturbations which get frozen on super- 
horizon scales, entropic perturbations could decay if the mass of entropic fields is higher than 
the Hubble scale. Especially, the moment when the inflaton path turns (e CT ^ 0) could come 
much later than the moment when entropic perturbations are generated, and it is possible 
that by the time e a ^ 0, all the entropic perturbations have already decayed away. If that 
happens, we will get N S Q S 0. 

In the context of a random potential, some of the entropic fields might be quite massive, while 
others are light compared to the Hubble scale. Only those light entropic fields contribute to 
the adiabatic mode and the number of such fields is D. It is entirely natural to have D S> 1 
in the context of string compactification. It is also possible that the mass of entropic fields 
fluctuates along the inflaton path. In such cases, if the mass does not remain larger than the 
Hubble scale for more than one e-fold, entropic perturbations may still survive. 



3. Inflation with a Random Multi-Field Potential 

One expects that the inflaton will follow a curved and random background trajectory in 
the Z)-dimensional field space. According to Eq.( 2.16j ), every time the inflaton makes a 



turn (e CT ^ 0), the perturbations along the e CT direction source SN. To determine SN, the 
evaluation of the integral in Eq. fl2.16 ) is crucial. In principle, such an integration can be done 
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given the knowledge of the scalar field potential and the inflaton trajectory in field space, and 
the result will be quite detail dependent. 

In this paper, we take an alternative point of view. We consider multi-field inflation with a 
random potential, and evaluate the power spectrum of the curvature perturbations generated 
by D canonical scalar fields. We focus on the scenario where D is large (D could be 100 
or larger); the scalar field interactions are complicated and can be strongly coupled, so the 
potential looks random (i.e., aperiodic). A typical scenario of such a scalar potential is the 
string landscape. In this scenario, the inflaton executes a Brownian-like motion in the D- 
dimensional field space even at the classical level. This type of random walk is different from 
the picture of stochastic inflation or eternal inflation, in which the random walk of the inflaton 
is driven by quantum fluctuations overwhelming the classical motion. Here the inflaton is 
mobile in the landscape. The random walk is due to classical scatterings in the random 
scalar potential and the resulting percolation. In our simplified scenario, the random classical 
scattering alone does not lead to density perturbation, which is sourced only by quantum 
fluctuation. However, the classical randomness of the potential does affect the mixing of 
adiabatic and entropic modes and thus will have impacts on the density perturbation by the 
end of inflation. In the next section, we shall focus on the scenario where the randomness 
comes from classical scatterings inside the random potential. The density perturbation comes 
from the quantum fluctuation on this classical background "random" inflaton trajectory. We 
will find that, when classical scatterings happen many times within one e-fold, the resulting 
power spectrum is relatively detail independent. In principle, a full analysis of quantum 
diffusion plus classical percolation will be interesting. However, since Sp/p ~ 10~ 5 is so small, 
such an analysis may not be necessary; that is, a perturbative treatment of the quantum effect 
on a classical background percolation is sufficient. 

It is not hard to imagine the following scenario. Start with a smooth multi-field inflaton 
potential suitable for slow roll inflation. Now, let us slowly turn on many tiny localized 
bumps with different heights and sizes randomly distributed in that potential. (One may 
include other features such as small dips, steps, valleys and ridges too.) In the limit of zero 
height, we get back the original potential for slow roll. The question is what is the effect of 
the bumps on the power spectrum. If the inflaton encounters only one bump during its travel 
in one e-fold, then one expects a feature to show up in the power spectrum. The size and 
shape of the feature in the power spectrum obviously depends on the properties of the bump. 
This has been studied in the literature already. Note that, order of magnitude wise, there 
are about 400 partial wave modes in the power spectrum per e-fold of inflation. Suppose 
that there is one data point per 20 partial wave modes. So, if the inflaton encounters more 
than 20 scatterings within one e-fold, we expect the distinct features due to the scatterings to 
smooth out in the data so one recovers a smooth observable power spectrum. It is in this limit 
that the statistical analysis presented here is fully applicable. As the number of scatterings 
decreases, our analysis becomes less reliable, but it also means that the likelihood of seeing 
some features in the power spectrum increases. We shall estimate the variance in the power 
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spectrum due to such scatterings. 



How likely is such a potential in the cosmic landscape? As explained in the introduction, a 
potential with small bumps is actually quite likely, since the barriers in the axionic directions 
are expected to be exponentially small. Of course, a pure axionic potential is periodic, 
not random. However, its coupling to other moduli (and axions) will lead to additional 
interactions that interfere with the exact periodicities, like impurities in a crystal. In fact, we 
hope that a detailed measurement of the power spectrum will tell us a lot about the particular 
patch of the landscape our universe has traveled over during the inflationary epoch. 

The randomness of the potential can come in different scales. One can approximate the travel 
of the inflaton as a long term deterministic drift plus frequent stochastic scatterings on short 
time scales. At the coarse-grained level, the inflaton will try to drift down the potential, 
dictated by the slope of the potential over large field scales. At the same time, the small 
scale randomness causes the inflaton to bounce around many times within one efold, leading 
to fluctuations in a on short time scales. We have in mind that these fluctuations are spaced 
by AiV e <C 1, i.e. the scatterings of inflation happen quite frequently within one e-fold. 

The coarse-grained deterministic motion down the potential results in a gradual change in 
the Hubble constant, and a small e parameter, defined as 



where H denotes the time average over the fluctuations within one e-fold. As the inflaton 
moves down the potential, e changes slowly, dictated by the rj parameter defined in the usual 
fashion by rj = e/(He). 

In the case with no or infrequent scatterings, we have e = a 2 /(2H 2 M 2 l ). By analo gy, we can 
define the coarse-grained speed (a) at which the inflaton drifts down the potential by 



In the following discussion, we shall always differentiate the drift speed {&) from the Brownian- 
like motion speed a. Here, a is the speed at which the inflaton travels between scatterings. 
We shall discuss two scenarios in the following sections. First, the weakly random scenario, 
where the drift term dominates over the scatterings. The inflaton undergoes frequent and 
mild scatterings, and after each scattering, e CT changes direction slightly. We expect {&) ~ a. 
Second, we also have in mind the strongly random scenario where the scatterings change the 
direction of e CT a lot. In this scenario, the stochastic scatterings dominate over the drift, and 
we generically expect (a) <C a. 

To distinguish the weakly and strongly random scenario, we define a parameter Xi 



e = 



2M 2 H 2 



=> (&)=HV2eM pl . 



(3.2) 



X 




(3.3) 
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X ~ 1 for the weakly random scenario and x ^ 1 f° r the strongly random scenario. We note 
that (a) measures the large distance properties of the potential while x ( or o") measures the 
short distance properties of the multi-field potential. 

A common feature for both the weakly and strongly random scenarios is that entropic per- 
turbation modes continuously feed into the adiabatic mode, resulting in the super-horizon 
evolution of the comoving curvature perturbations £• 

The field space distance traveled by the inflaton in the down-hill motion may be estimated 
by relating e to the overall tilt of the potential on large field scale V^U. Approximately we 
have 



2 



v 



2 



Using (a) = \ple HM p \ yields a downhill field length of 



N e V2~e Mpi « ^^-N e M^ . (3.4) 



However, on top of this deterministic down-hill motion, the inflaton also undergoes Brownian- 
like motion in the D dimensional field space due to scatterings on short time scales. 

To get a concrete idea, we can model the random walk using the Langevin Equation 



dt " 3H 



+m = v(<p)+m , (3.5) 



where all the vectors are defined in the field space. In general, \v\ > (a). To simplify the 
model, we have identified the drift velocity 

and £ is a random vector satisfying 

(C 1 ) = , (3.6) 
(£ 7 (ti)e J (t2)) ^2D IJ {$)5{t 1 -t 2 ) . (3.7) 

with D IJ ((j)) is the diffusion tensor. Here we have assumed the white noise property for £. 
This assumption is good if the scatterings take place at different sites in the landscape which 
presumably have uncorrelated properties and the typical time scale of the scattering itself is 
much less than the characteristic time between the scatterings. 

We can further work out D IJ . We assume the characteristic time between scatterings is At, 
and the typical field space distance between two scattering sites is A<p. Upon each scattering, 
the change of the inflaton position in the field space is given by 

8$=?At, (3.8) 
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and 5(f> can be decomposed along the direction of v and orthogonal to that, 

5$ = 5$ ± + 5$\\ , (3.9) 
If we assume that the typical scattering angle is 0, in the weakly random case with <C 1, 

\'2 



D-l ^ x ' V At At 

2o2 



Pi 2 At 

S^Sih-h), (3.10) 



D- 1 

<tyj[(Wf (fa)) * (A^) 2 ((l - cos 6) 2 ) eie J v s(j±-^ 

~ (A</>) 2 4 At Jfa - i 2 ) . (3.11) 

Here we define d 1 ^ = 6 — e^e^, with e v the unit vector along the direction v. We see that 
in the weakly random case, the longitudinal fluctuation is negligible. 

Since the drift velocity is also changing directions, we need to average over all directions of 
v, this will change S 1 ^ into 5 IJ (D — 1)/D. If we further use the relation Eq.(|3.8|) and the 
definition of D I J in Eq.(p.7|), we get 

£ /J -^fV J , e«i. (3.i2) 

If we introduce the diffusion constant A defined as D IJ = X5 IJ , we have 

e«:. 

In the strong random scenario, when randomly samples the interval [0, it], we have 

{54>l{t x )5<f>i{t 2 )) ~ { ^ )2 _f (sin 2 Q)5{ J 6(t 1 -t 2 ) , (3.14) 
(<tyj[(*0<tyjf (t 2 )) - (A^») 2 At(cos 2 6) 6(h - t 2 ) . (3.15) 

Note that the longitudinal fluctuation depends on cos(0) instead of (1 — cos©). The reason 
is that when 9 C 1, A(/>cos0 ~ Acj), and the fluctuation is given by A<^>cos0 — A0 = 
(1 — cos0)A</>. However, when takes all possible values from to it, A(/)cos0 itself is the 
longitudinal fluctuation. 

Averaging over all directions of v, we have 5{ J 5 IJ {D - 1)/D and -> 5 IJ /D. The 

diffusion tensor is 

D IJ - { ^/ J [(sin 2 0) + (cos 2 0)] * . (3.16) 
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In summary we have the diffusion constant 



e«,, 

A ~ , 9 ~ 1 . (3.18) 

2DAt ' v ; 



Equivalently, one may consider the time evolution of the probability density distribution 
P(4>, t) of the inflaton position in the multi-dimensional field space, which is described by the 
Fokker-Planck equation, 

dP 
~dt 



-V- 



P 



+ d!dj 



D IJ {4>)P 



(3.19) 



In the simple case, v is almost a constant vector and D IJ = A 5 IJ where A is the diffusion 
constant. In the strong random case, with G random over [0,7r], we identify 

\ _ (A0) 2 
AD At ' 

With these simplifications, the Fokker-Planck equation can be easily solved. 

Starting at <j> = at t = 0, the probability distribution of the inflaton at time t = N e /H is 
given by 

, —* . /A , , d ( \0-vt\ 2 \ / ,N e Y% ( \<f-vt\ 2 H\ 
P( M = (4,A f ) t exp (-^L J = (*rA-j exp J 

where the second expression is the probability distribution after N e e-folds of inflation. The 
amount of wandering explored by the inflaton is given by 

(<?) = 2DA§ 

As N e increases, the inflaton explores a larger field space. Here the field space explored in N e 
number of e-folds is roughly 

D + l 

(&)t(^)^~(&)(DX)^(^.j 2 . (3.20) 

Inflation ends when the inflaton reaches a big slope, or tunnels to a large dip in the potential. 
Another example is when the inflaton reaches a cliff (followed by tachyon condensation as in 
hybrid inflation or brane inflation) . Knowing the structure of the landscape (or some region) 
will go a long way to see if sufficient number of e-folds is natural. We know that the field 
distance explored by the inflaton goes like y/Ng, so naively many e-folds may be quite natural. 

We see that that the properties of the potential is captured in e, D, A and x- I n terms of 
independent parameters, we can trade A for At. We should point out that the treatment 
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of diffusion and drift in this section is very preliminary. In a realistic landscape, the drift 
velocity may not be a constant. The scattering angles many not follow a uniform distribution 
and the distribution may even vary over different parts of the landscape. The typical step size 
A(j) should also come with certain distributions and may be different in different directions. 
All these subtleties require a more careful treatment in the future. 



4. The Power Spectrum 

As we have discussed in the previous section, the small scale randomness of the scalar potential 
causes the inflaton to undergo Brownian-like motion due to frequent scatterings within each 
e-fold. As a result, the background inflaton trajectory makes random turns. To get an 
intuition on the behavior of the comoving curvature perturbation £ in the presence of these 
random turns, we can first look at Eq.(2.8). On the right hand side of Eq.(|2.8D, |e CT | is now a 



random variable. Every time e ff / 0, ( gets a random jump. So the D-dimensional Brownian- 
like motion in the field space translates into the one dimensional random walk of the scalar 
quantity £ through Eq.( |2.8[) . In fact, one can visualize this equation as the Langevin Equation 



for the one dimensional random walk [21] 



Using the 5N formalism, the power spectrum can be evaluated as 

P( = ^(NjQ'NjQ- 1 ) = ^ ((N a Q a N a Q a ) + (N S Q S N S ,Q S )) (4.1) 

In the above equation, the summation over all the s directions is implicit. If we further use 

H 2 



(QvQs) = , (QaQa) ~ 0/) , { , 



H 



2 



(N S Q S N S ,Q S ) = (N S N S ,)(Q S Q S ) = (N S N S ,)^6 SS , , (4.2) 



then we get 



P C = ^2 (N* + (N S N S ,)5 SS ') . (4.3) 

The appearance of (N S N S /) in the power spectrum is a particular feature of our scenario, and 
deserves more clarification. As we have seen earlier, 

fi>E 9/7 

N s = dt — eielj. (4.4) 
Jt* a 

When we calculate the power spectrum Pg, we take the variance over different Hubble patches. 
Due to the initial quantum fluctuations, different Hubble patch follows slightly different paths, 
and they make turns of slightly different angles in the field space. We emphasize that the 
difference in e CT among different Hubble patches, is sourced by quantum fluctuations Q s frozen 
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at i*. In this sense, we can talk about (N s N s r), even if N s appears to be a classical quantity. 
Note that N s is always accompanied by Q s in the form {N S Q S N S >Q S '), and we have factored 
out the quantum fluctuation (Q S Q S >). Such an approximation amounts to assuming that N s 
and Q s are weakly correlated, which is probably true in our Brownian like motion. Since Q s 
is frozen at the horizon crossing time t*, and N s received contributions every time inflaton 
scatters after t m , we expect that after a few scatterings, the correlation between Q s and 
N s becomes very weak. However, it was exactly because of the entropic perturbations Q s , 
that different Hubble patches will have slightly different N s , and (N S N S >) will be part of the 
resulting power spectrum P^. 

As we have explained earlier, e CT / only when the inflaton makes a turn in the field space. 
We assume the characteristic frequency of such turns is AiV e , i.e. there are l/AN e turns per 
c-fold. We then introduce 

(ei( tl )e J a (t 2 )) = JL&it^/Sih - t 2 ) , (4.5) 

where Sj~j = 5jj — ejej. 2 is the variance of the angle turned in the field space, it is the 
variance among different Hubble patches when they make slightly different turns at roughly 
the same time. Typically we expect < 1, since if different Hubble patches behaves too 
differently at each turn (large 0), 5N will be large by definition, and density perturbations 
can not remain small. 

Using the relation e*/ e*f = dj-j, we can write 

(N S N S ,5 SS ') = ^ _ N: s , (4.6) 



t* 



In the above definition of N* s , x is the time average of x over one e-fold, it appears because we 
integrate over the history of inflation, x it self may have fluctuations due to the Brownian-like 
motion, however x should be weakly time dependent. N* s denotes the effective number of 
e- folds from t* to ig. In the limit when both i9 and x are slowly varying, iV* ff approaches the 
actual number of e-folds iV*. The large scale observations of CMB temperature fluctuations 
today corresponds to 50 < iV* < 60. 



The power spectrum now reads 

P c (k) - 



H 2 H 2 

^Tx + ^ K *. 



(4.8) 

k=aH 



The first term comes from the pure adiabatic perturbation at horizon crossing. It depends on x 
which may exhibit fluctuations due to the Brownian-like motion. The second term containing 
$iV* ff ~ D0 2 N*ff/AN e characterizes the amplification of the power spectrum by entropic 
perturbations after horizon crossing. The factor DQ 2 describes the number of dimensions the 
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background random walk can turn to and the overall magnitude of the turning angle in field 
space. The factor N* s /AN e (^> 1) is the effective number of steps the random walk takes 
towards the end of inflation. The $N* S should be smooth since it depends on %■ The full 
shape of the power spectrum could exhibit some oscillations on top of a smooth background. 

We will first look at the power spectrum in the weakly random case. The inflaton undergoes 
mild but still frequent scatterings while moving down the potential. The stochastic scatterings 
do not cause sharp changes in e CT , and we expect •& and e to vary slowly too. Then we have 
X ~ 1, N* s ss N*. The power spectrum now reduces to 



k=aH 



The spectral index can be calculated using the standard method, 

d\nP c # , 

n s - 1 = - ; ps -2e - 77 , (4.10 

dlnfc ' l + N*ti V ; 

We see that if $N* > 1, n s - 1 = -2e- 77- l/N*. If -9N* <C 1, then n s - 1 = -2e- 77- 0. In 
either case, the correction to n s due to the multi-field effect is no greater than — l/N* ~ —0.02. 
We conclude that in the weakly random case, the power spectrum Eq. Q4.gp could easily 
resemble that of a single field slow roll model, and can easily agree with observational data. 

Second, we continue to look at the strongly random case, where both \ an d 7? vary with 
time. We expect that the long term drift speed (a) is much smaller than the Brownian-like 
motion speed &, i.e. % 3> 1. During the Brownian-like motion, & may fluctuate, leading to 
fluctuations in x, which will eventually impart oscillations on the pure adiabatic part of the 
power spectrum. Such oscillations in could be different from those generated by sharp 



features in the potential! 22, 23], in that they may appear incoherent with random amplitudes 
and phases. Since the background Brownian-like motion undergoes many steps within one e- 
fold, the fluctuations in should appear on all observable angular scales (with better chance 
of detection on small scales). The details of the oscillation in P^ cannot be fully captured 
by Eq. fl4,8| ) and numerical analyses are needed to calculate the evolution of £ around horizon 
crossing carefully]^]. However, the variance of the random oscillation might be estimated 
using Eq. Q4.8p . To do that, we first averaged out the fluctuations in P^ and introduce the 



averaged power spectrum defined as 

H 2 



8vr 2 ex 



P dk) = ^-j-:,[l+^K s } , , (4.11) 



k=aH 



The random oscillation about this average power spectrum can be characterized as 



y/HPt-ty) yCT (412) 

The variance of x can be directly related to the short distance properties of the multi-field 
potential, which is an interesting question to explore further. 
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The random oscillation of the power spectrum is a generic feature of our model, and may be 
falsifiable by future observations. If one looks at the un-binned TT CMB power spectrum 
one clearly sees that it fluctuates for about 10% (for 100 < / < 800). However, the error 
bars of each data point is also roughly 10%, so it is hard to claim any fluctuating features in 
the TT power spectrum based on current data. After binning a few tens of multiple moments 
into one data point, the current data analysis shows a relatively smooth TT power spectrum 
over all angular scales. On the other hand, from / ~ 100 to / ~ 1000, there are about 2 or 
3 efolds of inflation, so roughly there are 400 multiple moments per e-fold. If the inflaton 
encountered more than 20 scatterings per e-fold during inflation, then the current analysis 
based on the binned power spectrum will have smoothed out the fluctuations. In order to 
see such fluctuations in the power spectrum, it will be necessary to perform an analysis with 
the resolution of one or a few multiple moments. At current stage, such an analysis is still 
limited by statistics and systematics. With more data coming in the future, we expect to 
detect or falsify such primordial fluctuations in the power spectrum, if the error bars become 
substantially smaller than the primordial fluctuations. Furthermore, if such fluctuations are 
physical and appear in the TT power spectrum, the same fluctuations should appear in TE 
and EE power spectra as well, which will provide consistency checks on their existence. 

The evaluation of the spectral index is subtle in this case. As we are mainly concerned with 
the overall tilt of the power spectrum, not that imparted by the fluctuations in x, the spectral 
index should be derived from the averaged power spectrum P^. 

din A l-mN% „ 

""■Mi- "* s lk- (414) 

We have introduced new parameters rj& and r\ x to characterize the time dependence of i? and 
X- t]$ and r] x could be either positive or negative, unlike i9 which is always positive. 

Last, we want to mention that since the tensor perturbations are decoupled from the scalar 
perturbations, the tensor mode power spectrum is unchanged, while the scalar power spectrum 
gets modified. The tensor to scalar ratio r will be given by. 

( H\ 2 P T 16ex 
P T = 8 — , r = ^- = . (4.15) 

We see that the conversion of entropic modes into the adiabatic mode suppresses r by a factor 
of l/(i9iV* ff ). At the same time, r also gets enhanced by x ^ 1- Overall, r gets modified by 
a factor of xl (^-^Cff) compared to the single field slow-roll scenario. It is possible that this 
factor is greater than 1, e.g. one may have D ~ 100, AN e ~ 0.01, ~ 10 4 , N* s ~ 50, but x 
can be as large as 10 6 , if the course-grained potential is flat enough. The current observation 
bound on the tensor mode r < 0(0.1) gives, 

— ^—7 < 0.01 (4.16) 
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The analysis of the parameter space is more involved in the strongly random scenario. It is 
instructive to ignore rj x for a moment (assuming x varies slowly with time) and focus only on 
the last term in n s , which captures the effect of entropic perturbations. Below we will show 
two specific cases where some definite constraints can be drawn from observational data. 



1. \r)tfN* s \ » 1 and "&N* S » 1. The multifield effect on n s approaches 



i + m: s 



•& ~ \m\ • (4-17) 



For self-consistency, \rj^N* s \ 3> 1 requires \rj$\ S> 1/N* S . For r)$ < 0, i? decreases with 
time and N* s < N* » 60, leading to ^ <C -0.02. For < rfr < 1, we have N* s > N*. 
We generically do not expect N* s to be much different from N*, which ranges from 50 
to 60. This leads to 770 » 0.02. Recent measurement of CMB||] gives n s = 0.960^11, 
so we could easily accommodate —0.1 < 77,9 < 0.1. With a large tensor to scalar ratio 
r, the upper-bound in r]$ can further be relaxed. 



2. |r/^iV* ff | < 1 and $N* S » 1, then we have 

1 " ^eff 



# % ^ (4-18) 



In this case, since r]$ <C 1, varies slowly and we expect iV* ff ~ A^*. We may expect 
some running of the spectral index introduced by the —1/N* S term in n s , 



*a- ~ - A (".is) 



d In fe Af cff 



The WMAP+ACBAR dataQ has shown that dn s /dlnk = -0.037+0^23 is consistent 



with current observations, which implies an lower bound A^g > 10. We thus require 
that the scalar potential on large field scale do not get tilted too much, so that e does 
not increase too fast towards the end of inflation to give N* s <C 10. On the other hand, 
we also require entropic perturbations to dominate over adiabatic perturbations in the 
power spectrum to introduce such a scale dependence. If both conditions are met, our 
scenario could lead to detectable running of the spectral index. 



In this section, we have seen that the background Brownian-like motion can impart two 
features on the power spectrum. The super-horizon evolution of £ leads to a term proportional 
to $N* S in the power spectrum, which characterizes the conversion from entropic modes to the 
adiabatic mode. This contribution is continuous because the inflaton turns quite frequently 
in the field space within one e-fold. A second effect is that the small scale randomness of the 
potential leads to a "refraction index" x, which mainly affects the power spectrum during 
the epoch of horizon crossing. We have seen that in the weakly random regime (x ~ 1), the 
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power spectrum resembles that of the single field slow-roll case, with up to 2% more red tilt. 
While in the strong random regime (% 3> 1), the phenomenology could be quite rich. The 
tensor to scalar ratio can be enhanced and a large negative running of the spectral index may 
be observable. At the same time, the fluctuations in x could impart (high frequency) random 
oscillations on the pure adiabatic part of the power spectrum. The details of such oscillations 
might not be picked out by the CMB data. However, the variance can be estimated based on 
the short distance properties of the scalar potential. 



5. Non-Gaussianity 

Roughly speaking, possible sources of non-Gaussianity can come from three epochs: during 
the inflationary stage, during the onset of (p)reheating and after the universe has entered 
into the hot big bang radiation era. In the following discussion, we will first consider the 
generation of non-Gaussianity during inflation using the 5N formalism. This non-Gaussianity 
turns out to be small because the integration over the inflationary epoch tends to wash out the 
fluctuating contributions, and the growth of £ after horizon crossing suppresses the non-linear 
parameter /nl after proper normalization. 

Second, we will look at the mechanism to generate non-Gaussianity at the onset of (p)reheating 



in our scenario, which was first studied in Ref.[26, 27, |2£]]. This mechanism can take place 
in, for example, hybrid inflation, where tachyon instability ends inflation. The appearance of 
tachyons towards the end of inflation is quite natural in stringy inflationary scenarios. The 



abrupt ending of inflation is important for this scenario. Compared to Ref.[26, 27, \2^, we 
may gain a factor of D 2 or Z) 3 due to the large number of fields participating, provided their 
effects act coherently. 

Although the best known mechanism to generate non-Gaussianity in the radiation epoch, 
namely, the curvaton model|2£j, can be easily incorporated into our scenario, it is not a 
specific prediction of our model. So we will not discuss this possibility here. 

5.1 Non-Gaussianity Generated During Inflation 

The 5N formalism allows us to calculate the cubic order correlation of £ in a relatively easier 



way [30, 31 1. The nonlinear parameter /nl is defined as 

(C(k 1 )C(k 2 )C(k 3 )) = (2^) 7 ^|| (^/nl) + k 2 + k 3 ) . (5.1) 

Note that we normalize /nl in terms of the averaged power spectrum P^. 
If we write 

C = N I Q I + \Nu (Q V - <Q V» + . . . , (5.2) 
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following Ref.|31|, we can calculate the three-point function as 
(C(k 1 )C(k 2 )C(k 3 )) = (N I N J N K )(Q I (k 1 )Q J (k 2 )Q K (k 3 )) 



(5.3) 



+ ^(N I N J N KL )(Q I (k 1 )Q J (k 2 )[Q K Q L ](k 3 ))+ perms... . (5.4) 

The term in Q5.3| ) contains the bispectrum, and the term (5.4) is the product to two power 
spectra. The bispectrum can be parametrized as 



(Q / (k 1 )Q J (k 2 )Q A '(k 3 )) = (2^) 7 5(J» j0 
and we can calculate A IJK using the standard method 



Ekf (H 



2tt 



A 



UK 



(5.5) 



(Q / (k 1 )Q J (k 2 )Q^(k 3 )) 



dt < [Q I (k 1 )Q J (k 2 )Q K (k 3 ),H int (t)} 



(5.6) 



Here Hi nt is the cubic interaction Hamiltonian given by|L5|, |33|] 



int 



J d 3 x a 3 ^ 5jj (Q I Q J Q K + 2Q J d 2 Q J d- 2 Q K ) 



+ 



1, 



j^Vjj + -V IJK Q I Q J Q K 



(5.7) 
(5.8) 



The interaction terms involving Vjj and Vjjk are usually ignored in the single field slow-roll 
scenario, as they are higher order in slow-roll parameters. However, in the multifield case 
such assumptions might break down. In particular, Vijk could be large along the entropic 
directions. Furthermore, there exists small scale randomness that can cause Vjj to become 
momentarily large (3> H 2 ) during the background Brownian-like motion. Therefore, we need 
to include both the Vjj and Vijk interactions when calculating the bispectrum. 

Following the standard method, the three point correlation function due to the interaction 
terms (15.71) is, 



A. 



UK 



perm 



II s: J K 

dr V-rr^ e iKr [kjkl(l - ik 3 r) + 2k 2 k 2 {\ - ik 2 r)} + c.c. 



16H 



(5.9) 



The interaction terms (^J^) gives 



UK 



A 1J , 



perm 



dr 



?KllU 



+ 



32H 48H 2 



J — — JJ(1 - ihr) + c.c. 



(5.10) 



In the weakly random case, the above integrals can be estimated as 

■ b I 5 JK 



A 



UK 



16H 



auk 0KVU Vijk 
€r >' ~ 32 H ^ A8H 2 



(5.11) 
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where we have assumed that <fr, r]jj and Vijk are all slowly varying and we were not explicit 
about the momentum dependence. 

We point out that in the strongly random case, the careful evaluation of Eq.( |5.9| ) and Eq. fl5.10| ) 
is essential. Since rjjj and Vijk could fluctuate due to the background Brownian-like motion, 
there might be resonant enhancements in Ar e L [jl^]- However, by the nature of randomness, 
the fluctuations in r\jj and Vijk must come at a range of frequencies with random amplitudes 
and relative phases, and the effects may not add up constructively and may just self average 
out. If it happens that there are a few frequencies that dominate the resonance, A.*i K might 
be enhanced. We leave such a possibility open and more detailed analysis is needed. 

Given the power spectrum and bispectrum of Q 1 , the SN formalism tells us how to reprocess 
{Q 3 } to get (C 3 ). Unlike the single field case, where £ remains constant once leaving the 
horizon, here £ will evolve even on super-horizon scales, and the evolution is exactly dictated 
by the e-fold derivatives Nj and Njj in the 5N formalism. In order to estimate /nl> we first 
need to know these derivatives precisely. 



From Section EO, we know that 



w , = _£ e? _ r d(r ^M (5 . 12 ) 

a J a * a da 



We now continue to calculate Njj. Taking one more derivative in Eq.( 5.12| ) gives 

Njj = -dj (TW - + —^(dj«*) 

\a J a a da 



2H 
a 



-L d *v« d^- 2 !. d ^ B -{^) • (513) 



The expressions for Nj and Njj simplify a lot if we rotate it into the and e s directions. If 
we define N a = Nje^, N a p = Nije^ei, with a and representing either a or s, we have 

N a = ~- = —±= (5.14) 



a 



N s = -"r d aJ^^e- (5.15) 

VWt V ex da 3 K ' 

»- - ~ s ° (f ) = ^ ■ (5 ' 16) 

N - - -* (!) - vfetN • (5 ' 17) 

W ,„ = » r^M^^e-J + fd.^« , (5.18) 
V2e*x* J a * V ex da da J a « a e 



where we have used Eq.( [A.12 ) in deriving N ss i above. 
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We can now estimate the magnitude of non-Gaussianity. We first look at the bispectrum 
term (|5.3[), which leads to 



= ~3 • (5 - 19) 

We first consider the weakly random case, the A J e JK term contributes 

/NL ~ 2MN~Nj5") ~ "12(1 + ^) ~ "°-° 3r • (5 - 20) 
This term turns out to be proportional to the tensor to scalar ratio r, and using the constraint 
from the tensor mode Eq.( f4.16| ), we see that the above /nl is negligible. 

The term A e i K contributes 

, 5 {vijN 1 ^) 5 V IJK {NjNjNk) 

JNL ~ o /Ar ,r FT7A9 + 



8 (iV/iVj(5 /J ) 2 12 H 2 {NtNjS 13 ) 2 



r^> 




(5.21) 



In the above evaluation, we have taken N*^ 3> 1 and kept only the leading order terms. We 
have also assumed that the background Brownian-like motion gives (N a ) = and {N s N s iN s h) = 
0. Therefore, among all the Vhk's, only V acra and V ass i are left. We have defined new pa- 
rameters 

„ - f J *h y - v ^ 6 ± 

— D 1 ' ass — jj i 
We note that even if rj ss and V ass could be locally large in the multifield scenario, their 
contributions to /nl are still suppressed by l/(i?iV* ff ). Especially the tensor mode constraint 
Eq.(|4l6D further requires that exVss/($N* s ) < 0.01r/ ss . We generically do not expect large 
/nl coming from Eq.( |5.21 ), 



Again we should emphasis that the above results Eq.( |5.20|) and Eq.( 5.21 ) only apply to the 



weakly random case, where the scattering of inflaton is mild, although still being frequent. 
We still have entropic modes continuously feeding into the adiabatic mode. However, this 
conversion effect generically suppresses /nlj as such an effect enhances (£ 2 ) and (C 3 ) equally 
by a factor of #-/V* ff . Therefore, /nl ~ (C 3 )/(C 2 ) is generically suppressed by l/(i?iV* ff ). In 
the strongly random case, however, possible resonant effects might enhance /nlj which is not 
captured in our results here. 



Second, we look at the term Q5.4|) arises from the second order expansion of 5N and it involves 
the product of two spectra, and gives 

5 NjNjN KL 5 IK 5 JL 



/] 



NL 



6 (iY/7Vj5 /J ) 2 
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r, + Vx x 2 , 2^ x , 4^ + 4i?| s ^ 



fj ss is the time average of r] ss over one e-fold. We have also defined new parameters 



e(t) X (t) ■&* J ' J t . e(t) X (t) ff i 

Similar to iV* ff , N* and iVJ are just weighed averaged number of e-folds. In the limit that e, 
X, i? and f/ ss are constants, we have iV* ff = N£ = N%. 

Unlike /nl from the bispectrum, our result in Eq.( 5,22| ) applies also to the strongly random 
case. We now look at how large this /nl could be. We assume N* s ~ ~ iV|. More general 
cases can be analyzed similarly and the result remains the same qualitatively. We find that 
in Eq.( |5.22| ), we can at most get one term that contributes /nl ~ V-^cff' an< ^ an °ther term 
that contributes /nl ~ fjss/^- We will discuss them separately below, 



1. The 1/N* S term. As we have seen in the discussion on the power spectrum, scale 
invariance, i.e. the constraint on dn s /dln/c, requires N* s > 10, so this term contributes 
/nl < 1. 

2. The fj ss /•& term. In order for the large scale entropic perturbations not to decay away, 
we need the mass of entropic fields to be well below the Hubble scale, which implies 
fjss *C 1. One may expect that a small i? can enhance /nl; however, ■& is also bounded 
from below. In order for the entropic perturbations to have an effect, we need $N* S S> 1, 
which leads to i? S> 1/N* S . In the optimal case, we can take N* s ~ 50, fj ss ~ 0.1, then 
i? S> 0.01. We may have a chance to get /nl ~ 1 if we choose •& ~ 0.1. However this is 
still too small to be detected by future experiments. 



In this section, we have estimated the non-Gaussianity in our scenario. We conclude that /nl 
will generically be smaller than 0(1) in the weakly random case. The conversion of entropic 
modes into the adiabatic mode generically suppresses /nl by a factor of l/(#iV* ff ). Large 
/nl could possibly arise in two situations. First, in the strongly random regime, /nl might 
be enhanced resonantly. Second, one may change the properties on the background random 
walk, e.g. by introducing the correlation of e CT at different times or at higher order, then one 
could expect to have /nl 3> 1- We leave such possibilities for future study. 



5.2 Non-Gaussianity Generated During the Onset of (P)reheating 

Non-Gaussianity can be generated in hybrid inflationary models where inflation ends when a 
tachyonic field appears towards the end of inflationary epoch, which is quite natural in some 
scenarios inspired by string theory. In the string landscape with multi-fields, this can also be 
likely. Of course, both a and the tachyon field couple to other fields. The moment we are 
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interested in is when the tachyon field first appears, when a = oe- Let us focus on 2 specific 
terms in the SN expansion, 

5N = --Q° + -d s a E Q s + In ss , (q s Q s ' - (Q S Q S ')) +... (5.23) 

Let ifs be the entropic field with positive mass (squared) and couple to a, either directly or 
via the tachyon T. Let U(a, T, Lp s ) be the potential. Since 5Lp s will decay outside the horizon, 
due to their masses, m 2 ~ d 2 U/dLp 2 s (in the diagonal basis), we have 

Slp s ~ ^k s , k s ~ e -^ N °l™ 2 . 

where k s is the decay factor of the perturbation 5Lp s . If we further assumes that d s OE is 
dominated by the moment when tachyon appears, we have 



, p s = d s a E 

The non-Gaussianity generated is estimated in Ref.[^7]]. The Lp s fields can be non-Gaussian 
due to their interactions and this non-Gaussianity can be transferred to £ via the linear term 



in 5N (15T23|) . 



(( 3 )-^(^5Lp s ) 3 )\ tE (5.24) 
This yields, according to Ref.fl27f, 

/NL ~ — j= — rjn ^2 (.Ps K s){Ps>K>s>){Ps"Ks")Us S >s" , (5.25) 

where e* and H are evaluated as horizon crossing. Since ~ 1, we expect it to be much 
bigger than e*. With typical k s <C 1, we expect very small /nl generated here. Although the 
sum involves D{D+ V){D+ 2) /6 ~ D 3 terms, generically they may cancel among themselves to 
a large extent. It is also likely that the lightest Lp s (with the largest k s ) domintaes. Although 
we do not expect a large non-Gaussianity in general, the U ss ' s " terms may act coherently 
to make a big contribution to /nl- This is entirely possible since the contribution does not 
involve a time integration which tends to wash out random or fluctuating effects. 

The non-linear term also contributes to the non-Gaussianity. Because the non-linear term in 
5N is local in spacetime, they will yield a contribution of the local form 

/nl ~ ^4/2- E(^^)(^'4)C^| tj5 (5-26) 

^E s>s' 

This term can go like D 2 . Although we do not expect this local term to be too large generically, 
a large local non-Gaussianity is entirely possible. 
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6. Discussions 



So far in this paper, we have considered a random potential where scatterings in the potential 
lead to inflaton motion resembling that of a random walk. We analyze the impact of this 
classical percolation property on the density perturbation sourced by quantum fluctuations. 
In a more realistic situation, we expect quantum diffusion effects on the motion of the inflaton 
as well. By this, we mean fast tunnelling, resonance scatterings as well as quantum hopping. 
Here, quantum hopping is hopping over a low barrier, due to the presence of the Gibbons- 
Hawking temperature T = H/2n, so one may treat this as a thermal effect. In this scenario, 
we assume that the inflaton is not trapped at any classically stable site that will lead to 
eternal inflation. Rather, the inflaton is mobile in the random potential and can move freely, 
via classical percolation or quantum diffusion. The argument for this property is based on the 
multi-field feature and discussed in Ref.||. There, the multi-field feature (D 3> 1) is crucial 
for the mobility of the inflaton. Here we see that the resulting magnitude of the density 
perturbation requires many scatterings during a single e-fold of inflation, consistent with the 
overall picture that emerges. 

Instead of the decay rate per unit volume T, it is convenient to introduce the tunneling 
probability in one Hubble time, the dimensionless parameter 7 given by 7 = T/H 4 , so 7 
is a measure of the nucleation rate relative to the expansion rate of the universe. We are 
interested in situations where 7 > 1. Let the probability of an arbitrary point remaining in a 
false de Sitter vacuum at time t be P(t). Then P(t) behaves as P(t) ~ exp(-^-jHt). Thus 
the lifetime of the field in a false vacuum is estimated as 

3 3H 3 , . 

t F ~ = (6.1 

b 4ttHj 4vrr V ' 

If 7 3> 1, then Htp <C 1. Let us assume that inflation is dominated by repeated fast tunneling 



events. Following Ref.[19, 20|, or in the 5N formalism, we can easily estimate the amplitude 
of density perturbation 

C = SN ~ H5t ~ Ht p (6.2) 

COBE normalization is £ ~ 10~ 5 and thus tp ~ 10 _5 i^ _1 which implies that there are 
roughly 10 5 steps during one e-fold. This is the key idea of chain inflation fllOfl . We do like to 
note that an axionic potential with an exponentially small barrier height and a small decay 
constant do have the right feature for such a scenario. 

Here, we assume that each site has a dominant decay channel for fast tunneling, so that the 
inflaton path is uniquely determined. Since 7 > 1, the phase transition completes rapidly 
(that is, the universe percolates). So the inflaton path mimics a classical path. This puts a 
tight constraint on the inflaton potential. Presumably, this condition can be relaxed some- 
what; that is, the intermediate steps can lead to different paths, though towards the end of 
inflation, the inflaton ends up in the same vacuum site in the potential. Such deviation in 
paths can lead to a new type of density perturbation, which must stay small to satisfy the 
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observational bound. It will be very interesting to study this density perturbation and what 
it tells us about the inflaton potential. 



Even with fast tunneling, a locally attractive site will trap the inflaton even for a little while. 
In high dimensions, it is harder to trap a wavefunction. So in some cases, the inflaton is not 
trapped by a local attractive site, but rather it scatters as a resonance over that site. In both 
cases, e or a will decrease and so £ will increase. It is also likely that £ fluctuates, with vales 
much bigger than its average W (£ 2 ) ~ 10~ 5 . As long as such fluctuations happen at a scale 
too small to be picked up by the CMB data, (say, 10 3 steps per e-fold), it will be difficult to 
detect this effect. 

Such a scattering (or fast tunneling) from a meta-stable site A to a nearby site B can release 
a substantial fraction of the vacuum energy to the kinetic energy of the inflaton. Since the 
next step happens rapidly, there is no time for the Hubble expansion or the inflaton decay to 
dilute its kinetic energy, so next step scattering or tunneling can go to a site C with vacuum 
energy density comparable to that of A and above that of B. This means repeated scatterings 
and/or fast tunnelings do not have to go downhill. This leads to the Brownian-like motion. 
However, over some number of steps, the Hubble damping will take place and the inflaton 
will move downwards. This leads to the drift velocity. So at the coarse-grained scale, the 
resulting picture can easily mimic the slow-roll scenario. 

Even though the lifetime is much shorter than the Hubble time, the motion of the scalar field 
does not really roll down along a smooth potential, but instead a bumpy potential. Usually 
the small bumps do not affect the amplitude of the power spectrum, but it may impart a large 
distinctive non-Gaussianity fli~2|| . In our case the tunneling/percolation time is much shorter 
than the Hubble time, which means that the "period" of the bumps is very short and the 
non-Gaussian features from bumps may be too fluctuating to be picked up observationally. 
However, by the very nature of randomness in the potential, it is entirely possible that a big 
bump (or dip) in the potential or a particularly strong scattering can cause an observable 
feature in the CMB. This might be the origin of the feature at I ~ 20 in the WMAP data. 
Such fluctuations may appear in the power spectrum as a function of the wave number (or I). 
They are imprinted primordially, and should not be confused with the statistical fluctuations 
in data analysis. 

One possible subtle issue we have not discussed in this paper is the effect of isocurvature 
perturbations after reheating. When the universe is reheated through the decay of a large 
number of scalar fields, the abundance of the created particles may depend on the initial 
value of the inflaton fields. The isocurvature perturbations can be characterized as the per- 
turbation on the ratio of the particle number densities between different species, e.g. between 
photons and dark matter particles. Recent observations from WMAP+BAO+SN have put 
a stringent bound on the amount of isocurvature perturbations allowed in the primordial 



power spectrum[24]. The isocurvature power spectrum cannot exceed 6.7% of its adiabatic 



counterpart. In our scenario, the power spectrum of adiabatic perturbation at the end of 
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inflation is proportional to $N* S ~ (D — 1)Q 2 N* S , while we also have entropic perturbations 
by themselves with total power spectrum proportional to (D — 1), so the ratio of entropy 
perturbations over adiabatic perturbations is l/(N* s Q 2 ) which is about a few percent with 
O < 1. The amount of isocurvature perturbations left after reheating could be well within 
the current observational bound. On the other hand, the reheating process may result in ad- 
ditional power in the curvature perturbation and possibly additional non-Gaussianity. This 
is certainly an interesting effect that deserves future study. 

7. Summary and Remarks 

Motivated by the possibility of inflation in the cosmic landscape, which may be approximated 
by a complicated potential, we study the density perturbations in the multi-field inflationary 
scenario with a random potential. Here the randomness can come in different scales. It is a 
result of the complicated interactions among the moduli/axion fields in flux compactification. 
The small scale randomness of the potential causes the inflaton to undergo a Brownian-like 
motion in the Z)-dimensional field space, allowing the entropic perturbations to continuously 
and randomly feed into the adiabatic perturbations. The large scale randomness allows the 
inflaton to slowly move down in the landscape and eventually end the inflationary epoch. 

To quantify the density perturbations in this scenario, we employ a stochastic approach to 
evaluate the two-point and three-point functions of primordial perturbations. We find that 
in the weakly random scenario where the stochastic scatterings are frequent but mild, the 
resulting power spectrum resembles that of the single field slow-roll case, with up to 2% 
more red tilt. The strongly random scenario, in which the coarse-grained motion of the 
inflaton is significantly slowed down by the scatterings, leads to rich phenomenologies. The 
power spectrum exhibits such primordial fluctuations on all angular scales. Such features 
may already be hiding in the error bars of observed CMB TT power spectrum and have 
been smoothed out by binning of data points. With more data coming in the future and 
improved systematics, we expect these features can be detected or falsified if one analyze 
the TT power spectrum based on each single multiple moments. Furthermore, detection of 
such primordial fluctuations in the TT power spectrum implies the same feature in the TE 
and the EE power spectrum as well. On the other hand the tensor power spectrum itself 
is free of such fluctuations and the tensor to scalar ratio is enhanced by the large ratio of 
the Brownian-like motion speed over the drift speed. In addition a large negative running 
of the power spectral index is possible. Non-Gaussianity is generically suppressed by the 
growth of adiabatic perturbations on super-horizon scales, and is negligible in the weakly 
random scenario. However, non-Gaussianity can possibly be enhanced by resonant effects in 
the strongly random scenario or arise from the entropic perturbations during the onset of 
(p)reheating if the background inflaton trajectory exhibits particular properties. 

The inclusion of quantum diffusion effect into this scenario remains to be worked out. The 
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best guess made above is that it will not change the main consequences of our analysis. We 
believe that inflation with a multi-field random potential can be systematically analyzed with 
solid predictions. The study in this paper provides a first step along this direction. 
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A. An Explicit Derivation of SN under Entropic Perturbations Q s 

We now trying to calculated the perturbation of e- folds due to the entropic perturbation Q s , 
by considering a short segment of the inflaton trajectory as shown in Fig.[l]. The trajectory 
given by the dashed line (CD) is the perturbed trajectory from the solid line (AB) under Q s . 
The two trajectories correspond to different number of e-folds due to their different length 
(curvature radius) and different a along the trajectory. The number of e-folds along the 
unperturbed trajectory is given by 

N = —Al , (A.l) 
& 

where Al = adt is the length of the trajectory in the field space. Under the entropic pertur- 
bation Q s , Al changes to 1 

Al-AeiQj, (A.2) 
where Ae a is the change of the vector e a along the path. 

We also need the change of a under Q s . Since on super-horizon scales, Q s perturbations do 
not change the energy density d s p = 0, and given that p = |cr 2 + V(4> 1 ), we have 

d s & = . (A.3) 
a 



On the other hand, using the background equation of motion Eq.(2.2) and Eq.( |2.3[ ) one can 
easily show that 

+ ^ e to = -^., (A.4) 

a z a 



The sign convention is that the perturbation along the direction e a corresponds to a shorter path. 
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which can be written in components as 

= e s ■ e CT . (A. 5) 

a 

So we have 

d s & = e s ■ e CT . (A.6) 

We can now evaluate <5A^ due to Q s , we have 

N S Q S = --AeiQj - ^{d s a)Q s Al 



^AeiQj-S-eiQ^l , (A.7) 



a a 2 



The path length Al can be expressed as 



Al = j^-\Ae a \ , (A. 



where cr/l^o-l is the curvature radius, and |Ae CT | is the angle spanned by Al. 
Combing Eq.( |A.7| ) and Eq.( |A.8| ), we get 



N S Q S = --Aei Qj - E^QjlpA . (A.9) 
a a |e CT | 

The rate of change in 5N is given by 

sn = ^aT = ■ (A - 10) 

Using the relation that £ = SN, we have derived Eq.([2.§|) from the 6N formalism, 

C = -—eiQj . 
a 

Our result here does not depend on slow-roll approximations. 

Another relation we want to derive is the variation of e CT along the entropic direction, which 
was used in deriving Eq.( |5,18l ). Starting from Eq.( A. 4 ), we have 



e? = + *W 
a a 



taking derivative along the e s direction, we get 

= ~ Y -f + V -f^i + ^(a) - ^leW, . (A.ll) 
Especially, if we project d s (ef) to the e s ' direction, we have 

ef, c> s (e?) = + ^d s (&) = -- [V ss , + (e s , • e CT )(e s • e CT )] . (A.12) 
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Figure 1: Perturbation of the infiaton trajectory under Q s . The solid line (AB) is the unperturbed 
trajectory, while the dashed line (CD) is the perturbed trajectory under the perturbation Q s . 

B. The Evolution of Q a and Q s During Inflation 



The evolution of the comoving curvature perturbation follows Eq.(2.8) 



a' a 



(B.l) 



We have switched to conformal time dr = dt/a. We further introduce a transition matrix to 
quantify the time derivative of the field space basis, define 

is antisymmetric due to (e^e m /)' = 0. Since we have already chosen e K along e' a , we get 



Now Eq. flB.ip becomes, 



e' a = @ aK e K , e aa = (a//t) 



a' 



(B.3) 



(B.4) 



The evolution equation of Q n can be written in a nice form if we introduce v a = aQ a and 
v s = aQ s . 



k 2 



v a -2 {@ aK v K )' - 2—@ aK v K = , z = 

z rt 



j 2 a \ t- > 

k - — ) 6 s i + n sl 



vi - 2Q sl v[ = , {s + a, I + a) 



(B.5) 
(B.6) 
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where the effective mass matrix /j, s i for v s is given by 

& = a 2 V sl - Q' sl + sk Q kl + A@lJ SK 5 lK . (B.7) 

The fact that v a couples only to v K can also be understood from Eq.( |B.4| ), which can be 
rewritten as 

C' = - 2 ^^ • (B.8) 

Since C = —v a /z, we have 

v' a = z'( + C'z = z'C + 20 aK v K 

z" z' 

v" a = z"U + z'W + 2 (& aK v K )' = —v a + 2-Q aK v K + 2 (Q aK v K )' 



If we separate the solution to Eq^BJj) into the homogeneous solution and inhomogeneous 
solution v a , Eq.( |B.8|) provides exactly v a on large scales (k 2 <ti aH). Since only v K appears 
in C, we expect that v a is sourced only by v K . The homogeneous solution v® corresponds to 
the initial value of 

In terms of the 6N expansion, 5N = N a Q a + N a Q s , v® provides Q a and integral effects of 
v K provides N S Q S . In this sense, we claim that Q a and Q s are uncorrelated here. Since the 
effective mass matrix fi 2 , is not diagonal, in principle all the Q s perturbations are correlated. 
However, if the background inflaton trajectory is random, fi 2 , is a random matrix, the ran- 
dom correlation between Q s will get washed out over time. So we claim that Q s are also 
uncorrelated among themselves. 
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